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I. INTRODUCTION 



In recent years a large progress in the mathematical theory of the BBGKY hierarchy for 
quantum many-particle systems is observed. A good example of such progress is the rigorous 
derivation of quantum kinetic equations that describe Bose condensate [lH3| . In the original 
works of Bogolyubov (4-6] the solution of the initial- value problem of the BBGKY hierarchy 
was constructed in the form of the iteration series. The same representation of the solution 
is also being used in modern works 7H9|. 

For the case of the Maxwell-Boltzmann statistics the solution was also constructed in 
the form of the expansion over particle clusters, which evolution is governed by the cumu- 
ants of groups of operators for the von Neumann equations [10] (or by reduced cumulants 



111]) . These solution ex pan sions were constructed on the base of the non-equilibrium grand 



canonical ensemble 
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12|. 



In this paper we propose the alternative method of the description of the evolution of 
quantum many-particle systems. States of such systems are described in terms of correlation 
operators, evolution of which is governed by the von Neumann hierarchy. Based on the 
solution of such hierarchy we define the s-particle (marginal) density operators and derive 
the BBGKY hierarchy [5j, that can describe the evolution of infinite-particle systems. 
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Using the solution of the von Neumann hierarchy for the correlation operators we derive 
the formula of the solution of the Cauchy problem of the BBGKY hierarchy in the form of 
the expansion over particle clusters, which evolution is governed by the cumulants of group 
of operators of finitely many Fermi or Bose particles. 

Usual iteration series representation of the solution can be obtained from the constructed 
solution by the use of analogues of the Duhamel formulas for some classes of interaction 
potentials. 

Let us outline the structure of paper. In section 2 we define the basic notions of the Fermi 
and Bose many-particle systems and introduce the evolution equations for the correlation 
operators. In section 3 we introduce the s-particle density operators based on the solution 
of the von Neumann hierarchy for correlation operators and derive the BBGKY hierarchy. 
In section 4 we construct the solution of the initial- value problem of the BBGKY hierarchy 
in case of the initial data satisfying the chaos property. 



II. THE EVOLUTION OF CORRELATIONS OF FERMI AND BOSE MANY- 
PARTICLE SYSTEMS 

We consider a quantum system of a non- fixed (i.e., arbitrary but finite) number of iden- 
tical (spinless) particles with unit mass m = 1 in the space M 1 ", v > 1 (in the terminology 



of statistical mechanics it is known as nonequilibrium grand canonical ensemble 12j), that 
obey Fermi-Dirac or Bose-Einstein statistics. 

States of the system belong to the space £ 1 (J 7 ^) = ©^ -2 1 ('H^ : ) of sequences / = 
(fo, fi, ■ ■ ■ , fn, ■ • •) of trace-class operators f n = f n (l, ...,n) 6 ^{^-t) and f e C, sat- 
isfying the following symmetry condition: / n (l, . . . , n) = f n (ii, . . . ,i n ), if . . . , i n ) G 
(l,...,n), where H is a Hilbert space associated with a sing le particle, U± = W n is the 
symmetric (or antisymmetric) tensor product of n Hilbert spaces H; J 7 ^ = ©^L Wn i s t ne 
Fock space over the Hilbert space "H. Spaces £ 1 (J 7 ^) are equipped with the trace norm 

oo oo 

ll^ls 1 ^) = Sll^Ufi 1 ^) = Tri,..., n |/ n (l,...,n)|. 

n=0 n=0 

We denote by £q t ne everywhere dense set in £ 1 (J 7 ^) of finite sequences of degenerate 



operators 



13] with infinitely differentiable kernels with compact supports. Note that the 
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space contains sequences of operators more general than those determining the states 

of systems of particles. Hereafter we assume that ~H = L 2 (M. U ). 

The Bose-Einstein and Fermi-Dirac statistics endows the state operators with an ad- 
ditional symmetry properties. We illustrate them on kernels of operators Let 
f n (qi, . . . ,q n ; q[, . . . , q' n ) is the kernel of the operator f n G £ 1 (L^ ,± ). In case of the Bose- 
Einstein statistics the kernel is a function that is symmetric with respect to permutations 
in each group of arguments: 

fn(0.1i ?2j ■•■) Qn] Qli H2i —i In) = /n(?7r(l) > ?tt(2) > •■•> 1ir{n)] %'(\) > 9tt'( 2 ) ' '"' ^r'(n))' 

and in case of Fermi-Dirac statistics the corresponding kernel is antisymmetric: 

fn(qi, ?2j — > Qn'i q[, ?2; —J ?n) = (-l) (|7r| + |7r 'vn(?7r(l)> Qn(2), — > ?7r(n)i > ^'(2) ) — J 9^r'(n))> 

where 7r G (5 n and n' G (5 n are the permutation functions, 6 n is a symmetric group, i.e. 
the group of all the permutations of the set (1, 2, n), | vr | = 0, = if the permutation 
is even and \ir\ = 1, \tt'\ = 1 if it is odd. 

We define the permutation operator p n : 2}{'H < ^ n ) — > £}(H® n ) in terms of kernels of 
operators by the following formula: 

(Pnfn)(qi, 92) — , q n ] q[, <?2> ■■■■> Qn) = fn(t, 9tt(1)) Qn(2), — , qn{n)'i 9l> ?2) — > In)' 

Let us introduce the symmetrization operator S£ '■ 2}(1-L® n ) — > and the antisym- 

metrization operator S~ : 2}(H® n ) — >■ 2}(H~) by the following formulas: 

7re6n 

The Hamiltonian of the system H = @^L H n is a self-adjoint operator with the domain 

V(H) = {tp = ©Vn e J% | ip n e v(H n ) e nt,E\\ H n^ n \\ 2 < oo} C On the 

subspace of infinitely different iable symmetric (or antisymmetric) functions with compact 
supports ijj n G Lq (W n ) C L 2,± (IB^ n ) n-particle Hamiltonian H n acts according to the 
formula (Ho = 0) 

n n n 

H^ n = --J2^n+J2 ZJ $(fc) (?ii.--->ft fc )V'« > (2) 

8=1 fe=l ii<...<«fe=l 
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where is a fc-body interaction potential satisfying Kato conditions 13], h — 2%h is the 
Planck constant. 

We describe states of the system with Hamiltonian (T5]) by the sequence g(t) = (0, gi(t, 1), . . . , 
g n (t, 1, . . . , n), . . .) G £ 1 (J 7 ^) of correlation operators, which evolution is determined by the 
initial- value problem of the von Neumann hierarchy [141 ]: 

^g n (t,Y) = -tf n {Y)g n {t,Y)+ (3) 
+ E E •• E (-KF Zr \z 1 ,...,Z lPl ))S±l[g lXil (t,X l ), 

P:Y = \J i X i ZidXi, Z|p,cX| P |, XlCP 
|P|^1 Zi^0 ^|p|^0 

^n(t,i r )| t=o =Pn(0,l r ) > n>l, (4) 

where Y2p is the sum over all possible partitions of the set Y — (1, . . . ,n) into |P| nonempty 
mutually disjoint subsets X { C Y, cXj * s a sum over au subsets C Xj, for f„ G 
£q(H^) C V{M n ) C the von Neumann operator J\f n is defined by 



and 



In case of the Maxwell-Boltzmann statistics this hierarchy was studied in the work 

We note that the relation between correlation operators defined by hierarchy ©-(jlD and 
density operators D(t) G has the following form: 



9n{t, Y) — D n (t, Y) + (-1) |PM (|P|-1)!^ Hd^X,), (5) 

where the sequence D(t) is the solution of the Cauchy problem of the von Neumann equation 
111- 

Further we consider a more general notion, namely, the correlation operators of particle 
clusters that describe the correlations between clusters of particles. 

We introduce the following notations: Yp = ({Xi}, . . . , {X|p|}) is a set, which elements 
are |P| mutually disjoint subsets Xj C Y = (1, . . . , s) of the partition P : Y = |ji=i -^Q- Since 
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Yp = ({Xl}, . . . , {X|p|}), {Y} is the set that consists of one element Y = (1, . . . , s) of the 
partition P(|P| = 1). 

We define the declasterization mapping 6 by the following formula: 

e{Y P ) = y. 

Let us consider the set X c = ({Y}, s + 1, . . . , s + n). It holds: 

9(X C )=X = (!,..., s + n). 



The relation between correlation operators of particle clusters g(t) = (0, {^}), 



gi+ n (X c ), . . .) G and correlation operators of particles ^ is given by: 



9i +n {t,x c )= (- i ) |phi (i p i- i ) !5 s ± + n n e n m^ z v)- w 

P:X c =Ui X, X,CP P>: e(X i )=\J J Zy Z„CP' 

Using hierarchy ([3]) from formula (EI) we derive the generalized von Neumann hierarchy 
for the correlation operators of particle clusters 

j t 9i+n(t, X c ) = -M s+n (X)g 1+n (t, X c )+ (7) 

|p| 

+ S f +n E E •• E (-KF Zt \z 1 ,...,Z lPl ))l[g lXil (t,X l ), 

P:X c = \J i X i , Zice(Xi), Z| P | c 0(X| P |), XlCP 
|P| > 1 \Zi\>l I^|P|I>1 

where X = (1, . . . , s + n). 



III. THE DERIVATION OF THE BBGKY HIERARCHY FROM THE DYNAM- 
ICS OF CORRELATIONS 

Let us introduce the s-particle (marginal) density operators using the correlation opera- 
tors that satisfy hierarchy (J7J) by the following formula: 

X 

F s (t, Y):=J2- Tr s+1 ,..., s+n g 1+n (t, X c ), (8) 

n=0 

where Y — (1, . . . , s), X c — ({Y}, s + 1, . . . , s + n). Series ([8]) is convergent in £ 1 (J 7 ^) if 

gi + n{t)e&\H% n ). 
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We show that evolution of marginal density operators defined by expansion ([8]) is governed 



by the chain of equations introduced by Bogolyubov 



4]. The following derivation is given 



for the case of a two-body interaction potential, i.e. the terms J\f® t with I > 2 are equal to 
zero. 

Let us differentiate both sides of expansion (jSJ) over the time variable in the sense of 



pointwise convergence in £ 1 (J r w ) and use equality (J7J): 

oo 

;F s {t,Y) = ^-Tr s+1) ... jS+n ( -N s+n {X)g 1+n {t,X c )+ 



d_ 



n=0 



+ S f+n E EE (-M£!(h,i2))9\x,\(t,X 1 )g lX2l (t,X 2 )}, 

P:X c = Xl \jx 2 iiee(x 1 )i 2 ee(x 2 ) 

where X = (1, . . . , s + n) and the operators Sf +n are defined by ([T]). Taking into account 
the fact that 

K +n (x) = K(Y) + K(x \Y) + E M2(M 2 )), 

hCY i 2 CX\Y 

and that for g 1+n (t) G £ l 1 ('H s+n ) it holds 

Tr s+1> ... tS+n (-Af n (X \ Y))g 1+n (t,X c ) = 0, 
we rewrite the last equation in the following form: 

- F ^ t ' y ) = E^ Tr ^w + n((-A^(r) + E E (-^£K*ni2)))gi + n(t,x c )+ 

n=0 ' heY i 2 eX\Y 

+ S f + n E EE (-MSihMhxrlfrXjgwfrXz)). 

Since operator Af s (Y) does not depend on variables s + 1, . . . , s + n and according to the 
symmetry properties of operators gi +n (t), we obtain: 

d °° 1 

J t Fs(t,Y) = -K(Y)J2- Tr s+l,...,s + n9l+n{t,X c ) + 

n=0 

+ E ^Tr s+1 ,..., s+ „ (n E(--^2 (h s + l))gi +n {t, X c ) + 

n=0 ' i&Y 

+Sf +n E EE (-^inK^i2))g\x 1 \(t,X 1 )g lX2l (t,X 2 )y 

P:x c = x 1 \jx 2 ii^(x 1 )i 2 ee(x 2 ) 



d 
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According to definition (jSJ), first term on the right-hand side of the equation is equal to 
(—J\f s F s ). Using the symmetry property of the product g\ Xl \(t, Xi)g\ X2 \(ti X 2 ) which is the 
consequence of being under the trace sign, we make the following rearrangement in the last 
term: 

E EE (-^M))9\x l \(t,X 1 )g lX2l (t,X 3 ) = 

P:(X c ,s+n+l)=X 1 UX 2 ne0(Xi) i 2 e6(X 2 ) 

n+1 

= E E E EH^* 1 ' ^)M*> Z )92+n-k(t, (X c , S + 1 + U)\Z). 

k=l Z C(s + l,...,s + n + l), ^ Yi ^ Z 
\Z\ = k 

As the result we obtain: 

±F s (t,Y) = -M S (Y)F S (Y)+ 

+ E -TW,.,*+n+l ( E^^S^' S + l ))92+n{t, X C ,S + U+ 1) + 
n=0 ' ieY 

1 n+1 

-jCiE E E Y.(- N inSi,h))g k {t,z) x 



n + 

Z C (s + 1, .. . ,s + 1 + n) 
|Z| = fe 

xg 2+n - k (t, (X c , s + l + n)\Z)j . 

Then we get trace over s + 1 index out of the sum and use the symmetry property to rewrite 
the last term: 

j t F s (t,Y) = -K(Y)F s (t,Y)+ 

X 

+Tr s+ i ^(-A/2(^ S + !)) E — Tr S +2,...,s+n+l (^2+n(t, A c , S + H + 1) + 
iGY n=0 
n+1 j 

+«S s +„+i ^ ( fc _!)!("'+ i.fc)! ^' s + n + 2-k,...,s + n+l)x 

xg 2+n -k(t, {Y}, s + l,...,s + n-k + i)\ 

Now, to finish the derivation we need an auxiliary fact. In view of the symmetry of 
operators gi +n (t) with respect to s + 1, . . . ,s + l + n variables under the corresponding trace 
signs, it holds: 

g 1+n {t, {1, . . . , a + 1}, s + 2, . . . , s + n + 1) = g n+2 {t, {1, . . . , s}, s + 1, . . . , s + n + 1)+ 



S 



n+1 , 

+S± +n+1 ^ k _iy^'_ k + 1 y 9k{t, s + n + 2-k,...,s + n+l)x 

xg2+ n -k{t, {1, • • • , s}, s + 1, . . . , s + n - k + 1). 
Thus, we obtain the following equality: 

j t F s (t,Y) = -Af s (Y)F s {t,Y)+ 

oo ^ 

+Tr s+1 5^(-A/S(f, s + 1)) ^ -Tr a+2 ,..., a+n+ i0 1+n (t, {1, . . . , s + 1}, s + 2, . . . , a + n + 1), 

i& n=0 

and according to definition ([8]), we deduce finally: 

j t F s (t, Y) = -K(Y)F s (t, Y) + Tr s+1 £(-JV2(i, s + + 

This equality can be treated as hierarchy of evolution equations for the marginal density 
operators and was derived by Bogolyubov in work 4] from the von Neumann equation for a 
system of fixed number of particles. For the case of a n-body interaction potential we derive 
the BBGKY hierarchy in the similar way. Its initial-value problem has the following form 

j t F s (t,Y) = -Af s (Y)F s (t,Y) + (9) 

oo -. 

+ E^ Tr *+w+« £ (-^!Z l+n) (Z,X\Y))F s+n (t,X), 

lb. \ / 

n=l ZcY, 

F s (t) \ t=0 = F s (0), s > I. (10) 



IV. ON THE SOLUTION OF THE CAUCHY PROBLEM OF THE BBGKY HI- 
ERARCHY 

Let us construct the solution of the Cauchy problem of the BBGKY hierarchy (I9"j)- (TTU|) for 
one physically motivated example of initial data, namely the case of initial data satisfying 
the chaos property. A chaos property means that there are no correlations in the system at 
the initial moment of time (t=0): 

<7(0) = (<7i(MmO,0,-) ; (11) 
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which means in terms of marginal density operators dHJ that 

F 8 (0,Y)=g 1 (0,{Y}), (12) 

where Y — (1, . . . , s). 

Let us denote by 2li+ n (£, X c ) the (1 + n)-order, n > 1, cumulant of groups of operators 

G n (-t)fn = e-i tH »f n ei tH ", 

which is defined by formula: 

X 1+n {t,X e ) := (- 1 )' Phl (l P l " !)' II G\e(z k)l (-t,e(Z k )), (13) 

P:X c =U fc ^ fc Z k CP 

where H n is the Hamiltonian f n G £ 1 ('H^), 6 1 is the declasterization mapping defined in 
section 2 and JT C = ({V}, s + 1, . . . , s + n). 

In case of the Bose or Fermi system of particles a solution of the initial-value problem of 
the von Neumann hierarchy (CO) for initial data ( fill) has the following form 

g 1+n {t,X c ) = % 1+n {t,X c )Sf +n 11 9l {0,i), (14) 

where 2li +n (£, X c ) is given by (lT3"j) and operators Sf +n are defined by ([I]). According to 
formula (O), from definition (181) we obtain 



F s (t,Y) = ^Tr s+l _ s+n Vl 1+n (t,X c )Sf +n J] 9l (0,i). 

n=0 ' ieX c 

In view of equality (fT2"j) . i.e. 

</i(o,i)= n Fi(o,j), 

finally we obtain 

oo ^ s+ra 

F s (t,Y) = ^Tr s+1 _ s+n Vl 1+n (t,X c )St +n fl^CM, 

n=0 ' i=l 

which is the solution of the Cauchy problem of the BBGKY hierarchy (l9l)-( lT0l) for the Fermi 
or Bose many-particle system with initial data satisfying the chaos property. 

For arbitrary initial data F(0) G S^X^y) = ©r^=o OL n ^(Ji^ )> the following theorem 
holds. 
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Theorem 1 If F(0) G (.F^) = 0~ =o a n £ x {U±), then for a > e and i 6 R i/iere ezzsfe a 
unique solution of the initial value problem /fPj)-/ f70j) given by series 



oo _ 

F s (i, 1, . . . ,s) = — Tr s+ i i ... ia+n 2li +n (t,X c )F s+n (0,l,...,s + n). (15) 



ni 

n=0 



For initial data F(0) E £^ (J^) e ^(F^) expansion < T73)) zs a strong solution and for 
arbitrary initial data from the space £q(F^) it is a weak solution. 

In case of the Maxwell-Boltzmann statistics a solution of the BBGKY hierarchy (l9l)- (fT0|) 



for initial data from the space of trace-class operators was constructed in [10], |16|, [17] 



V. CONCLUSION 

In the paper the marginal density operators were defined by the means of a solution of 
the von Neumann hierarchy (J7J) for correlation operators by formula (JHJ). This definition 
allowed to construct the BBGKY hierarchy for the Fermi and Bose particles on the base of 
dynamics of correlations in the space of sequences of trace-class operators. 

We note that using the definition (jEJ) of marginal density operators it is possible to justify 
the BBGKY hierarchy in other Banach spaces that contain the states of infinite-particle 
systems, contrary to the definition of the marginal density operators in the framework of 
non-equilibrium grand canonical ensemble [ill ]. 



In the paper it was also defined the solution ffl5|) of the BBGKY hierarchy ([9]) for the 
Fermi and Bose particles with a n-body interaction potential. Such solution is represented 
in the form of the expansion ( IT5"|) over the clusters of particles, which evolution is governed 
by the cumulants ( Fl3l) of the groups of operators for the von Neumann equations. These 
cumulants for the 
Boltzmann case 



ermi and Bose particles have the same structure, as in the Maxwell- 
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